Electronic transport in an array of quasi-particles in the v — 5/2 non-abelian quantum 

Hall state 
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The Moore- Read Pfaffian v = 5/2 quantum Hall state is a p-wave super-conductor of composite 
fermions. Small deviations from v = 5/2 result in the formation of an array of vortices within 
this super-conductor, each supporting a Majorana zero mode near its core. Here we consider how 
tunneling between these cores is reflected in the electronic response to an electric field of non-zero 
wave vector q and frequency u>. We find a mechanism for dissipative transport at frequencies below 
the v = 5/2 gap, and calculate the q, u dependence of the dissipative conductivity. The contributions 
we find depend exponentially on \v — 5/2 1 1 . 



The v = 5/2 fractional quantum Hall state is ex- 
pected to be characterized by quasi-particles obeying 
non-abelian statistics. There are strong indications that 
this state is well described by the Moore-Read Pfaf- 
fian wavefunction which may be formulated within 
composite-fermion theory (each electron is bound to two 
flux quanta) as a p-wave superconductor of composite 
fermions (CFs) at zero magnetic field. Excitations in this 
superconductor are vortices carrying half a flux quantum 
and an electric charge of e/4, and fermions created in 
twos by breaking pairs with an appropriate energy gap 
0, The Bogoliubov-de-Gennes (BdG) equation de- 
scribing the fermionic excitations of a two-dimensional 
(2D) p-wave superconductor admits zero-energy solutions 
in the presence of well separated vortices, one solution 
near each vortex' core; These solutions are Majorana 
fermions 7, satisfying 7^ =7. As a consequence, the 
ground state is degenerate; For 2N well separated vor- 
tices, the ground state degeneracy is 2 N . The adiabatic 
interchange of two vortices induces a unitary transfor- 
mation within the subspace of degenerate ground states. 
Two such transformations do not necessarily commute; 
Hence vortex excitations obey non-Abelian statistics. A 
related spin model showing similar non-abelian excita- 
tions was recently studied by Kitaev Q. 

Experimental support to the Moore- Read theory is still 
needed. Relating the theory, and in particular the non- 
abelian nature of the quasi-particles, to measurable ob- 
servables, is a major theoretical challenge. Interference 
experiments may be a venue towards that goal 0, 0i • 

In this work we pursue a different method to probe the 
ground state degeneracy as well as some of the properties 
of the Majorana excitations, by considering the response 
of a quantum Hall system near filling factor v = 5/2 to 
an external electric field of wave-vector q and frequency 
w. In a fractional quantum Hall system at a filling factor 
v = 5/2 ± e (e <C 1), the density deviation from v = 5/2 
is accommodated by means of quasi-particles (vortices) 
whose density is 8en, where n is the density of electrons. 
For a perfectly clean system, these quasi-particles form a 
lattice, and when their density is large enough, tunneling 
between their cores should be taken into account. The 



degeneracy of the ground state is partially removed by 
this tunneling, and a band is formed with a width of 
the order of the tunneling strength. The tunneling also 
breaks the particle- hole symmetry of the localized 7i's. 

We study the electronic transport through that band 
for square and triangular lattices. We find that due to 
the existence of the band, there is a dissipative part to 
the conductivity below the v = 5/2 energy gap, with a 
unique q, ui dependence. This contribution to the con- 
ductivity, which does not involve a motion of the vor- 
tices, depends exponentially on lei" 1 / 2 , due to its origin 
in tunneling. There is a qualitative difference between 
the two lattice types. The square lattice is described by 
an effective massless Dirac Hamiltonian, while the trian- 
gular one shows a gap of a fraction of the band width. 
We calculate the dissipative part of the conductivity of 
the CFs using Kubo's formula ^ij, and then map it to 
the electronic conductivity by a Chern-Simon transfor- 
mation [i| (o- e ) _1 = (ct c/ ) _1 + f£e, (with e being the 
anti-symmetric tensor). For the square lattice, we find 
that the longitudinal and transverse CF conductivities 
are respectively 
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where rjn = oj — v^q . Here a is the lattice constant, Vq 
is the velocity characterizing the Dirac spectrum, and 1?, 
to be defined below, is related to the tunneling strength. 
For the triangular lattice, we find that the real part of 
the conductivity is 
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where ?? A = — 2 — ^-j-. As we explain below, the 
electronic conductivities are suppressed by a factor of to 2 
relative to the CF conductivities. 

There are four steps in the calculation leading to these 
response functions. First, we specify the Hamiltonian 
describing the array. This Hamiltonian turns out to be 
closely related to the Azbel-Hofstadter (A-H) Hamilto- 
nian HE!, describing electrons on a tight binding lattice 
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in a magnetic field. Second, we calculate the spectrum of 
the Hamiltonian. Third, we find how the system couples 
to gauge fields by expressing the density and current op- 
erators in terms of the Majorana operators 7j (with i the 
vortex index); we also present a physical picture of this 
coupling. Finally, we calculate the response functions. 

Based solely on the requirement of Hermiticity and on 
the relation 7.; =7], a lattice of well separated vortices 
is generally described by a tight binding Hamiltonian 



H = ^y%ij7i7j 
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where 7* are the Majorana operators satisfying {7^, jj} — 
Sij, and where i,j are nearest- neighbors lattice site in- 
dices. The tunneling strength t is real and positive. The 
matrix = ± is anti-symmetric and indicates the sign 
of the tunneling along the bond (i,j). While the freedom 
to redefine 7^ — > — 7^ makes the elements gauge depen- 
dent, the product of over bonds creating a closed path 
is gauge independent. We now show that this product is 
determined by a non-trivial phase a Majorana fermion ac- 
cumulates when encircling a plaquette, and give a simple 
formula for the effective flux per plaquette. This formula 
fixes the matrix up to a choice of gauge. 

In the absence of tunneling between vortex cores, the 
localized solution to a 2D p-wave BdG equation near a 
vortex embedded in a lattice of vortices is given by 

This is an approximate zero energy eigenstate of the first 
quantized 2D p-wave Hamiltonian i/edG ( see ) of 
an order parameter Ao(r) expif2(r; {R*}), where r is the 
2D-space coordinate and {Ri} are the vortices' positions 
and the phase f2(r; {Ri}) has the property of increasing 
by 27r around any closed path surrounding one vortex 
(clockwise). The phase appearing in the solution (@J is 
given by $ 4 (r; {Ri}) = Q(r; {R;}) + arg(r - Rj), where 
the first term originates from the order parameter and 
the second one originates from the p x + ipy pairing, which 
induces a relative particle-hole angular momentum. The 
point Pi is arbitrarily chosen close to the vortex core. 
The real wavefunction g(r) is localized at the vortex 
core. The tunneling matrix elements for nearest neigh- 
bors are purely imaginary, and are given by ±it where 



t = 
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tunneling strength, and where is the Hamiltonian 

in the absence of tunneling, of which Eq. Q is an exact 
zero energy eigenvector. For well separated vortices, t 
decreases exponentially with a ~ e -1 / 2 . 

To determine the matrix elements s^, we consider a 
Majorana operator hopping between n vortices along a 
closed path that forms a polygon whose edges connect the 
vortices. We show that there exists a non-trivial phase 
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related to this path, given by half the sum of the inte- 
rior angles of the polygon. The origin of this phase is 
in an interplay between the phase of the order parame- 
ter and the p-wave pairing. First we calculate the tun- 
neling matrix elements < Xil^BdclXj >— j ex P ^Pi j > 
where we use the tight-binding assumption to neglect 
the spatial dependence of the phase and explicitly set 
r = (R + Rj)/2 = Cij. For all bonds tij = t, while ipij 
is given by 

^ = \ X 3 vr!(i) ■ di + \ L 3 v arg(1 ~ Ri) ' dl 
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Varg(l-Rj) -dl (5) 



The first term depends only on the order parameter; 
it measures the change of the phase of the spinor due 
to vortices enclosed in the path. The second and 
third terms are the contributions to the phase due to 
the relative particle-hole angular momentum induced by 
the p x + ip y pairing; they measure changes in the di- 
rection of the path. Considering n tunneling events 
t n expi [V>ii» 2 + i>i 2 i3 + • • • + i>i n ii], the total phase is 
given by \ § VO(l) • d\+ 1 J27=i ^i, where Ai is the angle 
subtended by the path with respect to the i-th vortex, 
positive for anti-clockwise traversal. The first term gives 
a 7r-winding for each vortex enclosed in the path. For 
each of these enclosed vortices, Ai is given by minus the 
exterior angle, which can be written as — (27T — Jj), where 
Ii is the interior angle. For all other vortices Ai = h. We 
therefore get | / VO • dl + | J2i=i A i = \Hi I i, i.e. half 
the sum of interior angles of the polygon. This result is 
independent of whether the core of a vortex on the path 
is inside or outside of the polygon: if a path is deformed 
as to cross a core of a vortex, both the term related to the 
order parameter and the relevant angle Aj/2 acquire an 
extra 7r, and these two contributions cancel each other. 
For a general polygon of n vortices we get a phase of 
7rn/2 — 7r; Consequently, for a lattice whose plaquette is 
a polygon of n vortices we get n/4 — 1/2 flux quanta per 
plaquette. 

We note that for the A-H problem of tight-binding 
electrons on the same lattice with the same flux per pla- 
quette, the Hamiltonian is 



H 1 ' 
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The Hamiltonians © and JSJ share the same Harper's 
equation, their spectra are identical, but they differ con- 
siderably in the way they couple to gauge fields. Yet, 
there exist relations between their response functions. 

Our determination of the effective flux in © singles 
out then a chain of A-H type problems, one for each value 
of n, where the flux per plaquette is determined by the 
geometry of the lattice. There is a qualitative difference 
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between the triangular lattice, with an odd n, and the 
square lattice, with an even n; the former breaks time 
reversal symmetry in the effective A-H problem while 
the latter does not. The honeycomb lattice, for which 
n = 6, was considered in Q. 

In the next step we calculate the spectrum and eigen- 
vectors of the Hamiltonian After identifying the flux 
per plaquette, we choose a gauge which complies with 
it, commonly breaking translational symmetry. Trans- 
lational symmetry is restored by choosing a unit cell 
which contains an integer multiple of the flux quantum. 
The sites of the unit cell are numbered z = 1, . . . , s; In 
this way, we divide our lattice into s sublattices. We 
aim at finding an operator 1^ satisfying the equation 
[H, r^] = ET^ . We expand it in local site operators as 
= Ai7i ending with the following equation 



it ^ S%j \j — E \i 
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Using translational symmetry, the solution for each sub- 
lattice z can be written as A,; = e lk ' Ri A z (j) (k), where 
z(i) is the sublattice to which the site belongs. The 
equation for A z is given by 



fljEz'Aj 
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where the site i is an arbitrarily chosen lattice site that 
belongs to the z' sublattice. We denote the corresponding 
eigenvectors of H by \(°$ (k). This results in the following 
operators 



-.(a)t 



(k) 7i = ^AW(k) 



(9) 



which obey the usual fermionic anti-commutation rela- 



tions {4 a)t ,r^} = 6 aP 5(k~k>) and {1^,1^} = 
for positive energy modes. In terms of these opera- 
tors the Hamiltonian is diagonal and is given by H = 

Z^ka &kal k L k ' 

For the square lattice the A-H hamiltonian has half a 
quantum of flux per plaquette. We choose a gauge for 
which Sij = + along columns and has alternating signs 
between adjacent rows. Having translational invariance 
in doubled lattice vectors, we may split the lattice sites 
into four sublattices, numbered z = 1, . . . , 4. The Hamil- 
tonian H may be written in a 4 x 4 matrix notation as 

Hn = 2ta x (g> t z sm(ak x ) + 2ta x ® t x sm(ak y ) (10) 

In the limit \k\a — * the Hamiltonian 11U|) has a dou- 
bly degenerate gapless isotropic Dirac spectrum e^ a — 
sgn(a)uo|fc|, with a = +2, +1,-1,-2 and the character- 
istic velocity vq = 2at [ijj . The eigenvectors of Eq. itTUI) 
are given by 



(a)t r (/3)t! 



,(-2) 



,(2) 



lie' 



-i,l) 
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where e l9k = (k x +ik y )/\k\. 

For the triangular lattice the A-H Hamiltonian has a 
quarter of a flux quantum per plaquette. The Hamilto- 
nian in the sublattice representation is given by 



£T A = 2t ^ 6 sin(ai ■ k) + £ 3 cos(a 3 • k) 



(12) 



where £1 = 1(g) t x , ^ 2 — cr y ® r y , £3 = a y ® t z , and where 
ai = (ax — v / 3aj/)/2, &2 = iflx + y3ay)/2, a 3 = ax are 
the three lattice directions. There is a doubly degenerate 
spectrum indexed again by a, and the spectrum is [l^ 
£ka = sgn(a)e Aifc , where 



e A ,k = V2t\J 3 + cos(2ak x ) — 2 cos(ak x ) cos(V3ak y )(13) 

The spectrum is gapped, and there are two minima at 
ko = (±7r/3a, 0), around which the spectrum is quadratic 
£k +K,Q — sgn(a)v / 3i (l + ^a 2 K 2 ). The eigenvectors of 
Eq. (|12|l are given by 

1 



(1) _ ,(-2)* _ 



A,, = A 



= — (*B_ k ,-iS_ k) l,l) 

Jv k 



(14) 



A 



(2) 



A 



(-1)* 
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(-iB k ,-*B k> -l,l) (15) 



where is a normalization factor and 
Bk = eA(fc)/2^ + sin( ai -k) 
cos(a3 ■ k) + i sin(a2 • k) 

The coupling of the Majorana states of the Hamilto- 
nian © to an electric field is very different from that 
of the electrons in the A-H Hamiltonian ©, due to the 
particle-hole symmetry of the operators 7*. While each 
Majorana state (J2J is electrically neutral, when tunneling 
between vortex cores is switched on, a non-zero density 
of charge appears between the vortices. Projected to the 
subspace of the Majorana states, the density operator 
may be written as p(r) = Pij( r ) where 

p lj (r) = isijg(r - K t )g{r - Kjhaj (17) 

The operator i^cij has two eigenvalues ±1 which describe 
the sign of the charge mostly sitting at the center of the 
bond; However, the operators pij do not commute if they 
share a common Majorana operator, and consequently 
one cannot specify the charge at all bonds simultaneously. 
Furthermore, two nearest-neighbour spinors Xi an d Xj 
are exactly orthogonal < Xi\Xj >— due to a n phase 
difference between the overlap of the particles and that 
of the holes; however, they do support non-zero matrix 
elements of the charge operator < XiWzlXj >= isijfl, 
where $ = J g(r — ax)g(r). Consequently, the excitation 

rj, Q ^ carries a charge of -deka/t- 

Next we identify the current operator. At q = 0, the 
current is found using the identity 



j q=0 = i[H,d] = -m 



HjSjt 



R; — R?; 



7,7; (18) 
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where d = J r rp(r) is the total dipole operator. The sum 
2 j s ij s ji 7^ om y f° r sites i and I separated by a doubled 
lattice vector. The current may be transformed to k- 
space by inverting JHJ and substituting it into (|18|l . The 
q = current is a conserved quantity. To see that, we 
examine the commutator of j q =o with the Hamiltonian 



R, + R,: 



R, + R; 



[JJ,j q=0 ] OC ^ SijSjiSim \ — j ~;~„ 

ijlm 

which is described by paths composed of three consecu- 
tive bonds connecting the vortices given by I, m. All 
paths give zero contribution as they interfere destruc- 
tively with paths formed by starting from one of the ends 
and reversing the order of steps to the other end. 

At finite q we find the current by calculating p(q) = 
J r e lq r p(r) and using charge conservation. The current 
operator is, in momentum space, 



J(q) - 2^ ^'^Icq 1 k+q/2 1 k-q/2 



(19) 



where n 



a(3 
k,q 



ka/3 
(a)* 



^k+q/2 ' are *h e density matrix ele- 



ments of the associated A-H problem, and = ede^/t, 
Vk = <9k£fc are the charge and velocity of the quasi- 
particle respectively. For comparison, the longitudinal 
component of the current in the associated A-H Hamil- 
tonian is 



ka/3 



'k aW a)t c m 

K i VU c k+q/2 C k 



q/2 



(20) 



where for the relevant transitions J{ 1 



Jo 



»k-q/2j f or the square lattice and 
= ^|(3 + 3V3i) for the triangular lat- 



<|a (e^+q/ 2 - e 

t/i jh 

J l,-2 — J 2,-l 

tice near the bottom of the band. Over all, the matrix 
elements of the current operator in the Majorana prob- 
lem l|19|) are smaller by a factor of qa relative to those of 
the A-H problem (f2~0")l . 

Having calculated the spectrum and identified the rele- 
vant operators, the response functions of the array of Ma- 
jorana states is readily calculated employing the Kubo 
formula, with the results given by Eqs. 0J, (J2J- This 
result affirms the existence of dissipative conductivity, 
hence the flow of in-phase current, even at frequencies 
below the v = 5/2 energy gap. 

Two steps need to be taken to transform the composite 
fcrmions conductivities and (0) into the measurable 
electronic conductivity. First, the imaginary part of the 
CF conductivity, ip s e 2 /uj (with p s being the superfluid 
density of the CFs) , originating from the superconductiv- 
ity of the CF condensate, should be added to the calcu- 
lated real part. Second, the Chern-Simon transformation 
should be used to transform the CF conductivity into the 
electronic one. In the limit of small lu, these two steps 
result in 



Recr[j(q, oj) 



(— 

\2h Pt 



Reer 



of 



Reoj_(q, uj) 



\2hpt 



Re a 



cf 



(21) 



At finite temperature, assuming v$q <C u>, the conduc- 



tivity satisfies <r[T] = a[T = 0]sgn(o;) tanh 



ik B T- 



Before closing, we note that the same methods may 
be used to find the response functions of the associated 
A-H problems. For the square lattice the conductivity 



is Re (cr^j I, , cr 



8 h 



1/2 

1/2 I 

in 



The value of 



the conductivity at the q — > limit is a universal e /8H 
|l3| ; the origin of the universality lies in an exact cancel- 
lation of the dependence on vq due to the linear density 
of states oc uj/vq. The dependence of the conductivity 
on temperature is Op = | ^-sgn (uj) tanh 4 ^ T ■ For the 
triangular lattice the conductivity at the bottom of the 
band is again universal, 



3 el 

4 h 



In summary, we calculated the electronic response of 
an array of immobile quasi-particles of the v = 5/2 state 
to an electric field of non-zero q, lu, due to tunneling be- 
tween the Majorana modes at their cores. We found a 
contribution to the dissipative conductivity, Eq. I|21(l . 
that is of a unique q, u> dependence, and a strong expo- 
nential dependence on the deviation of v from 5/2. Our 
analysis neglected disorder, which will be discussed else- 
where. 
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